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Abstract 



We present for the first time a self-contained calculation of the QCD running coupling at fi- 
nite temperature and quark chemical potential, a a (T,n), based on a semiclassical background field 
method. The hard thermal/dense loop results on the Debye screening mass are recovered in a first 
approximation. The final result can be interpreted as the ordinary zero temperature running cou- 
pling, with momenta replaced by in-medium scales A»: at high density and zero temperature, the 
, quark scale is set by A* ~ 24.4 fx. At high temperature and moderate densities, the quark in- 

medium scale reads A 2 ~ [2.9f T] 2 + [1.91 fi] 2 , which resembles the naive phenomenological estimate 

' 1 Introduction 



A 2 = (^r) 2 + M 2 . 



First steps in the QCD phase diagram have been made possible recently by progress in lattice calculations 
at finite temperature T and baryochemical potential fis = 3/x, where fx is the quark chemical potential 
^ 121 El- Results on the equation of state of QCD are now available for small to moderate chemical 
potentials fi Si T c . Whereas the temperature region close to the phase transition is inherently non- 
perturbative and at the moment only treated within phenomenological models, both at fi = |SJ |H] 
and fi 7^ 0E], various perturbative resummation techniques apparently work quite well above roughly 
3 T c |UlllU| and allow an interpretation of the thermodynamics in terms of weakly interacting quasiparticle 
partons that carry the same quantum numbers as the fundamental quarks and gluons. An extension 
of these calculations to finite chemical potential is straightforward and will allow further tests of the 
models by comparing their predictions to the mentioned lattice data. 

However, since there is no scale around other than the renormalisation scale A in the chiral limit of QCD, 
all the rcsummations basically do is to fix the magnitude of the partition function. The variation with 
respect to temperature is solely caused by the temperature-dependent running coupling constant; the 
normalised, dimensionless pressure, e.g., reads p(T)/T 4 — f[ct s (T/ A)]. In all perturbative calculations, 
a s is assumed to take on the form of the common zero-temperature running coupling a s (k 2 ), derived 
from one- or two-loop renormalisation group equations. The momentum scale k is then usually chosen to 

*'Work supported in part by BMBF, GSI and by the European Commission under contract HPMT-CT-2001-00370. 
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resemble a typical thermal scale A*, e.g. the lowest non- vanishing gluon Matsubara frequency 2itT '. In 
[111 IT2] , we have derived that this is indeed a sensible choice for long wavelengths provided the different 
Matsubara frequencies for bosons and fcrmions are taken into account: 



a s , e ff(J , A) = 



a. (A) 
12tt 




(i) 



with 
and 



A^ = A g (2wT), where ,4 9 = exp(- 7 - 1/2 log(3/8) + 1/44) » 0.938, (2) 



A{ = ^/(7tT), where ^/ = exp(-7 + 1/2) w 0.926 (3) 

(stars * denote in- medium quantities). The separation of thermal in- medium scales for the quarks and 
gluons arises naturally in our calculation and is also expected phenomenologically on second thought, 
but has not been taken into account in actual perturbative calculations so far. 

In the presence of a quark chemical potential, it is not so easy to argue what the relevant scale might 
be. First, one expects that /i modifies, to one-loop order, only the quark part of the running coupling, 
hence A» = f(T,n). Second, the fermionic Matsubara frequency is not an intuitive guide line anymore 
since it becomes complex: u> n (T,n) — (2n + l)irT — ifj,. From |wo| 2 or, equivalently, identity l|4tj[l . one 
might argue that 

(A{) 2 = (7rT) 2 + M 2 (4) 

constitutes an appropriate in-medium scale for fcrmions. For \x T, the temperature would hence 
dominate the magnitude of A{ . On the other hand, one may define an average thermal momentum by 

d 3 k\k\f D (k- fi) 

(fe)* = ^ ? ~ 3.15 T + 3/i + 0(m 2 ) (5) 

d 3 kf D (k) 



for small /jl/T, and /ij(fc) = [exp(/c/T) + 1] . This result implies that T and /i are on an equal footing, 
in contrast to our first estimate above. 

Also note that the in- medium momentum (k)* for quarks will - on general grounds - always increase 
with the chemical potential. A larger (fc)„ (or, equivalently, A{), however, implies that the antiscreening 
property arising from the gluons will be even more compensated by the (ordinary) screening from the 
quarks than at fj, = 0, hence 

> 1 ^ 

a s , e s(J- ,0) 

The presence of a chemical potential might then invalidate perturbative calculations that work at zero 
li. Taken literally, for each fixed value of the temperature there exists a (very high) density where the 
screening part of the quarks will eventually overcome the gluonic antiscreening and reverse asymptotic 
freedom. 

It is useful to recall at this stage the hard thermal/dense loop (HTL/HDL) perturbative results on the 
effective charge. From the static limit of the longitudinal part Dx(fc°, k; T, /i) of the gluon polarisation 
tensor one finds (fc = |fc|) ^Hl 

a SjCff (fc 2 ; T, fi) = rr ffTT^f ,A = "57^ — 7 ( 7 ) 

llz,(0, A:; I , (i) m A D {T,n) 

fc 2 + fc2 
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with the Debye screening mass 



2 ,-r , I Nc , N f\ 2rr,2 , N f 2 2 / c s 

m D (T, n)=(— + -J-\gT + ^g fj, . (8) 



Note that eq.Q is both momentum- and medium-dependent. Since gluon and quark contributions to 
m D have the same sign, both contributions screen, and the ratio 

a"^(^;r,0) <1 f0rallfc ' (9) 
Correspondingly, nothing spectacular happens when T and/or \i are taken very large. 



Clearly, various physical arguments suggest different and partly contradictory statements about the 
nature of a s (T, /i), which justifies in itself a detailed analysis. Also, in view of the importance of a s (T, fi) 
for future calculations, it is desirable to have a more refined calculation that aids to set quantitatively 
the in-medium scale A{. In this paper, we will extend the formalism of refs.[nj[^], that has already 
proven to work very successfully at fj, = 0, to finite chemical potential. A great advantage of the method 
lies in the fact that the extension is unique and physically intuitive. Section 2 briefly reviews the basics 
of the calculation. In section 3, we treat the high density limit at T = 0, section 4 finally features the 
calculation of the physically relevant case of high temperature and fi SS T. 



2 Basics 

In |1 II I12j. we have extended the approach of refs.^3 E] to finite temperature. Instead of a loop 
expansion of the gluon self-energy, the thermal energy shift 

AE(H, T) = ~ [47r X (tf, T)] VH 2 - £(T) vac , (10) 

of the perturbed thermal vacuum to order a s , after applying an external static chromomagnetic field H, 
was evaluated at temperatures T 3> Aqcd ■ The calculation of AE involves a sum over all Landau levels 



u n ,k 3 ,s3 = \Jki + 2eH{n + l/2 + s 3 ), (11) 

weighted by the corresponding thermal occupation probabilities. Here, n labels the mode, &3 is the 3- 
component of the momentum and S3 is the ^-component of the spin of the particle under consideration. 
Identifying 2eH with the scale k 2 at which the system is probed, as at T = (e is the product of the 
strong coupling constant g times the charge number q that involves the structure constants of SU(iV c )), 
we have investigated the high temperature limit 

eH 

^2 « 1- d 2 ) 

Note that condition i|12|) is equivalent to the hard thermal loop (HTL) approximation in thermal per- 
turbation theory. Subsequently, we have extracted a temperature- and momentum-dependent dielectric 
permittivity e(fc,T) by use of the relation 

a s Mk, T) ee -J5L_ = ^J^, (13) 

which is also valid at zero temperature. Here, x{k, T) = \g + Xq 1S the sum of gluon and quark magnetic 
susceptibilities. An expression similar to l|lUfl . the pure quantum part, persists even at zero temperature 
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and has to be added to itl"U|) in order to obtain the total result, so Xtotai = Xo(k, A) + x(fc, T). The zero 
temperature susceptibility xo{k, A) has been calculated in |13U14| and results, using (|13fl . in the famous 
QCD running coupling expression: 



4^ X0 (H, A) = An [xg(ff, A) + x g(ff, A)] 



48tt 2 



HJVclog 



2eH 



2N f log 



2eH 
T 2 " 



(14) 



A corresponds to a zero temperature scale that characterises the "medium" and is ultimately identified 
with the renormalisation reference point of the coupling. We will be concerned with the in-medium 
part of the quarks, Xq(k,T, fi), in the following and add the vacuum contribution \o and the gluon 
contribution Xg(k,T) only in the very end. 



3 The high density limit 

We will first evaluate the high density limit, i.e. T — and 

eH 



< 1. 



For simplicity, we assume isospin symmetry in the following. The quark chemical potential acts only on 
the quark energy density, therefore we sum up all energy levels up to the Fermi surface: 

AE{H,fi) = & -El 

E f = ^ ^n,fc 3 ,s 3 [6(Q - fi) + 6(Q + fi)] and 

n,k 2 ,k 3 ,s 3 

u k [9(uJ k - fi) + 0(uj k + fi)] , 

ki,k 2 ,k 3 ,s 3 



F f - 



with 0(x) the Heaviside step function, and summation over flavours is implicit. Apparently, antiparticles 
do not contribute to AE. For the unperturbed case, uik = \k\ = k, and 

d 3 k 



E{ = V 



2N 



(27T)"- 



Ak\9(k-fM) 



N f , 



This expression has to be subtracted from 

E f V ^ 



4tt 2 



i i /o n •/ 



u- fi), 



(15) 



(16) 



s 3 =±l/2 n=Q 



where we have replaced the summation over the quantum numbers ki and k^ by an integral and a 
corresponding density of states, as outlined in JT]. Performing the spin sum, separating off the n = 
contribution and transforming to dimensionless variables, we find 



E f = V(eH) 



JJ_ 

4tt 2 



— Y / dx \/^ 2 + n + l d(\Jx 2 + - fi) 

77 n=o{ 



(17) 



This expression is still exact. To proceed, we have to make some approximations. As in we first 
simply replace the sum over n by an integral. Again, this step will turn out to be insufficient to capture 
all effects in e 2 , and we will give the complete result below, with corrections to the replacement taken 
into account. For now, we obtain after transforming to polar coordinates, 



E f = V(ftff) 



2 

^- + [dr vV + 1 0(Vr 2 + 1 - fi) 

VK 1 -K Z J 



(18) 
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Since we work in the high density limit, ft 3> 1, the integral is dominated by the large r limit. Approxi- 
mating \/ r 2 + 1 ~ r, wc find 



(19) 



The term c>c fi 4 is just Eq that has to be subtracted - an important check of the consistency of the result. 
Summing over flavours and charges, the final result reads 



47rx(H",M) 



9 2 N f m 2 
2eH 2tt 2 



2eH 



(20) 



with the Debye screening mass (JSJ at zero temperature and high density. With the help of ca. lfHfl) . wc 
recover the perturbative result Q, exactly as in ^J, when we neglect corrections to the density of states 
and the higher- lying Landau levels. 



We now take into account these corrections: first, we use the Euler-MacLaurin formula 
fin) = / f(n)dn + - [f(N) + /(0)] + - [f(N) - /'(0)] + ... 

n=0 { 

to get the density of states correctly up to order e 2 : 



(21) 



E* = V(eH)-^ + V(eH) 2 ^r 

V M7T 2 V ' 7T 2 



1 

24 



du- 



1 



2 J duu 2 \Ju 2 - 1 + i y 
1 

— /t) + ?i<5(w — /!)} 



Second, we evaluate the integrals, expand the results in the small parameter 1 /ft, 2 — 0(e) and keep all 
terms up to order e 2 . Then, 



Ef = V(eH)J^ + V(eH) 2 -2 



8(eHy 



8eH 48 48 



"2ei? 



(22) 



Despite the very different structure of the integrals, the final result bears a surprising resemblence to 
the calculation at high temperature and zero density: All terms of order \j? cancel, and the remainder 
shows the characteristic log structure. Hence, 



^ Xq {U^)^g 2 —L\ 0% 



V 2eH 



, B f = 2e 5/2 ~ 24.4 



(23) 



As in the high temperature calculation, the pre-factor of the log is exactly the same as in the zero- 
temperature case, so, including xo(H, A) of ea. ll4|) and trading 2eH for the momentum scale k 2 , the 
final result reads 



af (fc 2 ,A; M ) 



a s (A) 



1 + 



a a (A) 
12tt 



HJV c log 



A 2 



2N f log 



f[B f lA 



2 \ 1 



V A 2 



(24) 



What is somewhat surprising is the large number B> , which we expected to be C(l) from eq.lQJ. Note 
that the final result (|24|l is both k 2 - and /x 2 -dependent, so it defies a simple interpretation like "take the 
zero temperature, zero density running coupling a s (k 2 ) and simply replace k 2 by /i 2 ". The situation at 
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finite density is different from the high temperature case, where both gluons and quarks acquire a thermal 
average momentum, since it singles out the quarks. Furthermore, the large number B' suggests that the 
antiscreening of the gluons will be largely compensated due to the very large average momentum of the 
quarks. It is therefore not obvious when a s (k 2 ,^) will become so small such as to justify perturbative 
calculations. 



4 High temperature and moderate density 

Next, we explore the combined effects of temperature and chemical potential. The finite temperature 
gluon contribution 



[A?] 2 
2eH 



(25) 



has been dealt with in ^2 and is unaffected to this order by n, so we investigate the quark sector 
only. Again, we start from 

AE(H,T,fx) = E'-E*, 

E f = ^ Q n,k 3 ,s 3 + /£>(<*> + M)] and 

n,k 2 ,k 3 ,S3 

[/d(w* - n) +/c(wfe +//)] . 

h\ , &2 , &3 .S3 



The unperturbed energy is simply that of the non-interacting fermion gas, 

-7tt 2 1 M 2 1 M 4 ' 



E£(T,v)=VT 4 N f 
which is easily obtained from ea. ^58|) . since 



60 2 T 2 4vr 2 T 4 



E f Q = VT^Nf 



24 

lim F 5 (eAp/T)/e) 

TT Z e->0 



Now, in appropriate dimensionless units, 



(26) 



s 3 =±l/2n=0 



Ef(H,T, fJ ,) = VT 4 ^ £ / dx^x* + e 2 (n + 1/2 + s 3 ) f±(y/x 2 + e 2 (n + 1/2 + s 3 )), 



where 



2 2eH _ /I 

6 = mn and U = —, 



and we defined, for the sake of lucidity, 



fax) 



exp(.T — /x) + 1 exp(x + (i) + 1 



After performing the spin sum, 
E s = VT 4 



o 2-2 2 00 

(T e fi e — 



24 8tt 2 2vr 2 



n=0 



J] / c/aVz 2 + e 2 (n + l) fay/x 2 + e 2 (n + 1)) 



(27) 



(28) 



which is still exact. In the following, we explore the high temperature, moderate density limit, i.e. 

e 2 < 1 and /2 <, 1, (29) 
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and calculate an approximation to eci. (|28[) up to order e 4 /i 2 = g 2 (n/T) 2 . 

We proceed in close analogy to the high density case: first, we simply replace the sum by an integral. In 
our HTL/HDL approximation, we neglect corrections to the sum and the Landau levels to arrrive, after 
summing over the couplings, at 

F _ V h 2 \ ^ 52712 "fWX + F* 

E f-~2 VH I" 6" 2eH~^2eH\ +E °- (30) 

The quantity in curly brackets is already 4nx q (H, T, /i), which leads, when including the gluon sus- 
ceptibility 4irx g (H,T) in the corresponding approximation |1 lj and setting k 2 = 2eH, to an effective 
coupling 

a s>eS (k 2 ; T, fi) = 3 (31) 

k 2 

with the HTL/HDL Debye mass ©. Recovering again the perturbative result by exactly the same 
procedure as at high temperature or high density, respectively, is an important check of the internal 
consistency of the used method here. Our incorporation of the chemical potential in the calculation is 
therefore well founded, which adds credibility to the complete result of 0(g 2 ), to which we turn now. 

Following our previous approach, we first take into account corrections to the density of states according 
to eq.JHJ: 

E f = VT i ^ + ^ + ±[24F 5 (e,[l/e) + 2e 2 F 3 (e : fi/e)] + 

2 4 \ 

+ ±^ [2F 3 (e,/2/e)+e 2 ^(e,M/e)] - JL-H(e,fi/e)) , (32) 

using the expressions for the Fi and H, as defined in the appendix. Second, we now expand the 
subsequent integrals in the small parameter e with the help of eqs.lJSHJ, JS7J|, and up to order 
e 4 to arrive at the simple expression 

AEf(H,T,n) = AEf(H,T,»= 0) + VT 4 {-IM.^ , (33) 

with C(3) ~ 1.202... All terms of order e 2 /2 2 cancel identically. Taking into account Xo(^^) °f ea. (|14|l . 
we obtain the following form for the total quark susceptibility: 



4n Xq (T^)^g 2 ^L 



log 



f[A f 7TT} 2 \ , 7C(3)^ 2 



V A 2 J 2tt 2 T 2 



(34) 



which is H -independent, as in the high temperature, zero density case. This is a highly non-trivial 
result that only arises because of intricate cancellations of terms other than C(e 4 ). The persistence of 
scale-independence even in the presence of a chemical potential is a remarkable outcome. 



5 Conclusions 

We have presented for the first time a self-contained calculation of a medium-dependent QCD running 
coupling constant a s (T, /i) at finite quark chemical potential, based on the semi-classical background 
field method of refs.m^]. We have considered the two cases high density, zero temperature (i.e. 
[i 2 /H ^> 1 and T = 0) and high temperature, moderate density (i.e. T 2 /H ^> 1 and fi/T ^ 1), and 
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have calculated the induced energy shift AE? of the vacuum in the quark sector. In both cases, the 
corresponding HTL/HDL perturbative results were recovered by neglecting corrections to the density of 
states and the higher-lying Landau levels, as in the high temperature, zero density case. The approach is 
therefore internally consistent, and the inclusion of the quark chemical potential is free from ambiguities 
and sensible. 



Taking into account the mentioned corrections, all H- and hence scale-dependence drops out of the 
in-medium contributions once the vacuum (zero temperature, zero density) contributions are added. In 
addition, the running coupling becomes very simple in both cases, as at high temperature, zero density: 
it again follows from the zero temperature renormalisation group equations, with the momentum scales 
of the propagating particles replaced by appropriate in- medium scales A*. At high density, the gluon 
part remains unmodified, whereas the quark part has the form 

,2 2N f ,_([M] 2 \ 



■log 



a 48tt 2 \ A 2 
with the in-medium scale A{ = 2e 5 / 2 /x ~ 24.4 ji. 



At high temperature and moderate density, the final result that follows from l|34fl . including the gluon 
contribution xfi + Xg{T), can be written as 



-off 



1 



12tt 



llAUog 



\ A 2 



2N f log 



[Hi 

A 2 



(35) 



Here, the gluon in-medium scale is obviously the same as in the high temperature, zero density case, i.e. 

A^ ~ 0.938(2ttT). 

The quark in-medium scale receives a correction due to the chemical potential: 



Bf = 



which can, for fi/T ^ 1, be expanded as 

[A{] 2 ~ [2.91 T] 2 + [1.91 fi} 2 



7C(3) 
2tt 2 



0.43, 



(36) 



(37) 



This form of the quark in-medium scale A{ is surprisingly reminiscent of our initial estimate J3J . We have 
therefore been able for the first time to give a sound (though semiclassical) justification for approximate 
use of the naive phenomenological quark scale A 2 = (ttT) 2 + /x 2 in the zero temperature, zero density 
running coupling, although in our calculation the influence of the chemical potential is somewhat larger 
(which is, however, in accord with the result at finite density and zero temperature, see above). 



To get a quantitative feeling, we show in figure 1 the effective coupling g° s (T, /i) as a function of T/T c , 
for fi = 500 MeV. There, g s (T, /i) (solid line) is our exact result (|35|l . g s (T, 0) the same at /i = (dashed 
line). Apparently, the coupling strength increases, as anticipated, at finite chemical potential due to 
the larger average in-medium momentum of the quarks, and the effect becomes most pronounced at 
small temperatures. For comparison, gP hen (T, /i) corresponds to the phenomenological expectation, i.e. 
A? = 2ttT and A{ = y/(nT) 2 + (dash-dotted line). Applying the full result ijSEJ still leads to a 
coupling strength g s larger than one in the temperature region where lattice results on the equation 
of state at finite /x exist. It might nevertheless serve as quantitative input for actual perturbative 
calculations that will try to describe these data. 
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5 

T/T 
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Figure 1: The effective coupling gg S (T, fj,) as a function of T/T c 
MeV and n ~ 500 MeV. For explanations, see text. 



T c is taken to be 170 MeV, A QC d ^ 250 



6 Appendix 



To make this paper self-contained, we present here the formulas used to evaluate the e- and /i-dependent 
integrals, appearing, e.g., in eq. l|32[l . Since e <C 1 plays the role of a small mass term, the whole machinery 
of the small m/T expansion of thermodynamical integrals can be applied. Our approach correspondingly 
simply extends the calculations of |16| . where the expansions were carried out for bosonic integrals, to the 
fermionic case. We will be short on some details, hence, but comment on differences where appropriate. 
We start from 



9n(y,r) 



1 



r(n) 



dx x r ' 



1 



exp(y / x 2 + y 2 — ry) + 1 



(38) 



and 



fn(y,r) 



dx- 



r(n) J y/x 2 + y 2 



exp(i/a: 2 + y 2 - ry) + 1 



Here, Y{n) is the Gamma function. It is useful to define 



G n (y, r) = g n (y, r) - g n (y, -r) 



(39) 



(40) 



and 



F n(y, r) = f n (y, r) + f n (y, -r). 



(41) 



These integrals appear in thermodynamical expressions for an ideal gas of fermions with mass m at finite 
chemical potential. The corresponding dimensionless quantities are therefore y — m/T and r = /i/m. In 
the following, we will be concerned with the high-temperature limit of ea. H41l) . i.e. m/T <C 1 and /i < T, 
for odd n. In the bosonic case, the physical region is \r\ < 1, otherwise Bose-Einstein condensation 
phenomena have to be taken into account. For fermions, no such restriction exists, so r can be in 
principle arbitrarily large. The expansion of (|40|l (and, correspondingly, l)4ipl in small y will have the 
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general form 

oo , 

4 (y,r) =/„(», r) + X>*(r)(J) (I = F,G), (42) 

fc=0 

where /„ might be a non-polynomial function in y, and the a^(r) are themselves polynomials in r. If 
dfc(r) is of the form 

k 

a,k{r) = by, with b k ^ 0, 

i=0 

expression l|42|) also becomes an expansion in ry = p, = fi/T. To be able to capture all leading effects in 
we will therefore restrict ourselves in the following to the situation p, ^ 1. 



G n and F n obey the same recursion relations as their bosonic counterparts: 



With the initial conditions 



— Lrn+i — 

ay 


mF„+i 


, y 2 r 

H fn-l 

71 


y r 

Lrn-l, 

71 


d F 
ay 


r 

— G n -i 

71 


V F 

71 




(71 > 0), 


F n (0,r) 


2(1 -2 2 - 

71- 1 


-n) 

'C(n 1) (ti > 2 



(43) 
(44) 

(45) 



being the Riemann Zeta function, all necessary functions can therefore be derived from G\ and -Fx- 
In close analogy to ^Hl > we substitute the identity 



1 i °° 

1 1 ^ z r AR \ 

TI = 2~ ^ z 2 + (2ti + 1) 2 tt 2 ( ^ 



n— — oo 



+ (271 + 1) 2 7T 2 

in the integrands of G\ and F\ and integrate term by term. We start with G\ : 

oo 

G 1 (y,r) = -2j2S m (y,r), (47) 



m— 

where 



'■'» = 2 <<"' I dx i 2+ ^^L-Xf^^ V ( 48 ) 
[x z + y z (l - r z ) + uj^\ z + u>f n [2yry 



and the (dimensionless) Matsubara frequencies cu m — (2m + l)ir. Note that Q m has exactly the same 
functional form in the bosonic case, only there u) m = 2rmr. Next, the integrand in (|48|l is expanded in 
powers of y and integrated term by term. After that, the summation over m can be performed, yielding 
£ functions. The final result reads 



G 1 (7/,r) = -r 2 / + 2^^(-l) fc A[l-2-( 2fe + 1 )]C(2fc + l)(^ , (49) 



k=l 

1 

2 



where A\ — 1, A2 = 2r 2 + §, A3 = 3r 4 + ^-r 2 + These coefficients are the same as in the bosonic 



For F\, we need to introduce a convergence factor x e (0 < e < 1) in the integrand. In the end, the 
result is finite, and the limit e — > can be safely taken. We start from 

00 

F 1 (y,r)=T(y,r)-2j2^m(y,r). (50) 

m=0 
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Now, 



and 



f{Vi r ) = / dx 



= - - log 



(|)+0(e)' 



P m (y,r) = 2 dx x 



\J x 2 + y 2 e 

x 2 +y 2 (l-r 2 )+u 2 m 
[x 2 + y 2 (l - r 2 ) + u 2 m } 2 + u 2 m (2yrY 



The leading term in y is 



(51) 



(52) 



T&\y,r) = 2 / dx 



' x 2 + uj 2 n Lu^t e cos(7re/2) ' 



(53) 



Using f(l + e) = l/e + 7 + 0(e), 7 = 0.5772... being the Euler-Mascheroni constant, the summation over 
the Matsubara frequencies yields 



00 -. 
2 £^)( y ,r) = -i-7 + log(J), 



(54) 



m=0 



which removes the 1/e singularity of T{y, r). Higher order terms in y are straightforward, and the final 
result reads 



00 

F x {y, r) = - log (J) - 7 - 2^(-l) fc S fe [l - 2 -( 2fc+1 )]C(2fc + 1) 



2k 



(55) 



fe=i 



with the universal coefficients Bi = r 2 + |, B 2 = r 4 + 3r 2 + |, B 3 = r 6 + ^r 4 + ^ r 2 + 4- Here, 



Sfc oc r , therefore the first term in /x in F\ starts at order /j, y 



2„,o 



With cas. 1)43(1 . I)44|l and 1)45)1 . all other terms of interest can be worked out to any order in y. We quote 
the expansions that are of interest in this work: 



(56) 
(57) 





Fi(y,r 


) = 


-log 


F 3 (y,r) 


TT 2 

= I2 + 


4 


r 2 + 


/7tt 4 tt 2 


[r» - i" 


y 2 


3 


^120 + T 


2 


+ 16 





" 2 1" 




r + 2 


-5 + *(*)] 


7C(3) 




16tt 2 





y4+ ^C(3) 


' 2 1_ 




©; 


y 256tt 2 


r + 6. 


y 6 ) 



„ 2 3 
- 2r 2 - 7 + - - 
4 

(58) 

Note that a term like r 2 y 6 is actually of order (yit/T) 2 (m/T) 4 , so it is necessary to push the expansion 
to higher terms in y to capture all effects in m/T to a desired order. 

As in |TJ, for the evaluation of the derivative terms in ea . (|32l) we also need the small y expansion of 
integrals such as 



H(y,r) 



dx 



\j x 2 + y 2 1 + cxp(— \\J x 2 + y 2 — ry]) 



(59) 
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exp(-y/ x 2 + y 2 — ry) + 1 



+ (r — > -r) 



exp(-\/a; 2 + y 2 — ry) 
^exp(-\/a; 2 + y 2 — ry) + lj 



+ (r — > — r) 



(60) 
(61) 



For the evaluation of the second term in (|61|l we employ the same trick as in |TI] . Introduce a parameter 
a to write 



exp(ay/x 2 + y 2 - ry) 
ax ^^=^= h (r — > — r) 



(exp(ay x 2 + y 2 — ry) + I) 2 



H(y,r;a) = - 

Obviously, -ff (y, r; 1) is the sought quantity. Now we re-write 

1 1 



(62) 



H(y,r;a) 



d_ 

da 



dx 



^ x 2 + y 2 exp(a^/x 2 + y 2 — ry) + 1 



+ (r — > -r) 



Therefore 
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